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Abstract 



We report an apparent detection of the Cvi 4p to Is transition line at 459 cV, during a long-term 
enhancement (LTE) in the Suzaku north ecliptic pole (NEP) observation of 2005 September 2. The observed 
intensity of the line is comparable to that of the C VI 2p to Is line at 367 eV. This is strong evidence for 
' the charge-exchange process. In addition to the Cvi lines, emission lines from OVII, OVIII, NeX, and 

(•~^ \ MgXI lines showed clear enhancements. There are also features in the 750 to 900 cV range that could be 

ff^ ' due to some combination of FeXVII and XVIII L-lines, higher order transitions of O VIII (3p to Is and 6p 

0^ ' to Is), and a NelX line. From the correlation of the X-ray intensity with solar- wind flux on time scales of 

, about half a day. and from the short-term 10 minutes) variations of the X-ray intensity, these lines most 

likely arise from solar-wind heavy ions interacting with neutral material in the Earth's magnetosheath. A 
hard power-law component is also necessary to explain the LTE spectrum. The origin of this component is 
' not yet known. Our results indicate that solar activity can significantly contaminate Suzaku cosmic X-ray 

spectra below ~ 1 keV. Recommendations are provided for recognizing such contamination in observations 
Q ■ of extended sources. 

^ ' Key words: atomic processes — Sun: solar wind — Earth — interplanetary medium — X-rays: diffuse 

c/3 . background 

- ^ , 1. Introduction related with the solar-wind proton flux, which strongly 

^ ■ suggests the origin of the LTEs to be solar-wind charge 

H \ Snowden ct al. (1994) reported the existence of myste- exchange with Hi in the geocorona. Solar- wind charge 
rious X-ray contamination episodes in the i?OS'^T all sky exchange results in line emission from highly ionized 
survey data which they termed long-term enhancements ions (Krasnopolsky et al. 2004 and references therein). 
(LTEs). During an LTE, the X-ray counting rate in the ROSAT did not have enough spectral resolution to resolve 
lower energy bands as much as doubled on a time scale those lines. Spectral information on geocoronal solar- wind 
of 1-2 days. However, they could not find any correla- charge exchange was first obtained during a Chandra dark 
tion with other observational parameters, such as space- moon observation (Wargelin et al. 2004). The X-ray pho- 
craft position or look direction. New insight on LTEs tons detected in the direction of the dark moon are most 
was obtained by the discovery of X-ray emission from likely from this source. The emission spectrum could be 
comet Hyakutake (Lisse et al. 1996). Following the dis- described by a sum of Cvi, OviI, and OVIII K-lines, 
covcry. X-rays were detected from many comets (e.g., although the statistics and energy resolution were lim- 
Dennerl et al. 1997; Cravens 2002), and the emission ited. More recently, Snowden, Collier, Kuntz (2004) re- 
mechanism is now well understood as charge exchange ported time variation of soft X-ray intensity during the 
of solar-wind heavy ions with cometary neutrals (see XMM-Newton Hubble deep field north observation. The 
Krasnopolsky et al. 2004 for a review). Then Cox (1998) enhancement was correlated with solar- wind proton flux 
and Cravens (2000) suggested that solar- wind charge ex- variations. They detected C VI, O VII, O VIII, Ne IX, and 
change with neutrals in the Earth's geocorona and in the MgXI emission lines in the enhancement, 
heliosphere accounts for a part of soft X-ray background The importance of solar- wind charge- exchange emission 
below 1 keV. Robertson et al. (2001) showed that the is three-fold. Firstly, it will enable us to remotely study 
LTEs of the ROSAT all sky survey data were well cor- low density neutrals in geocorona, in outer atmosphere of 
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planets such as Jupiter, in interplanetary space, and in 
particular, around comets. Secondly, it can be used as 
a highly sensitive ion probe for the solar wind. Perhaps 
most importantly, it becomes a significant contaminating 
foreground in the study of the cosmic soft X-ray back- 
ground below ^ 1 keV. More than half of the soft X-ray 
background at these energies is considered to arise from 
hot gas in the disk and halo of our galaxy and in in- 
tergalactic space. Although such emission was detected 
early in the history of X-ray astronomy (e.g. Tanaka, 
Bleeker 1977), its origins and the physical state of the hot 
gas are not yet well understood. Geocoronal interactions 
with the solar wind produce at least a sporadic contamina- 
tion that must be avoided in observations of any extended 
source. Charge exchange with interstellar neutrals mov- 
ing through the heliosphere creates a more subtle contri- 
bution where much of the diagnostic rapid time variation 
is washed out by the large travel time of solar- wind events 
through interplanetary space, and the spectral lines arise 
from the same ions expected in hot interstellar plasmas. 
Lallement (2004) used a simplified model of interplane- 
tary charge exchange to estimate that essentially all of 
the minimum fiux observed near the Galactic plane in the 
1/4 keV (R12) band of the ROS AT sky survey could arise 
from this source (see also Pepino et al. 2004). We are far 
from an understanding of solar-wind charge exchange ad- 
equate to determine the true extent of this contribution. 

The X-ray Imaging Spectrometer (XIS) on board 
Suzaku (Mitsuda et al. 2006) has a significantly improved 
spectral line response function compared to previous X- 
ray CCD cameras, particularly below 1 keV (Koyama et 
al. 2006). Together with the large effective area of the 
X-ray telescopes (XRT; Serlcmitsos et al. 2006), this will 
open a new era for the study of soft X-ray background. In 
this paper, we report on Suzaku observations of a blank 
field in direction of the north ecliptic pole (NEP). We de- 
tected a significant enhancement of the soft X-ray flux 
lasting for ~ 10 hours. The enhancement is mostly ex- 
plained by increases in C VI through MgXI emission lines. 
During the enhancement, both C VI n = 2 to 1 and C VI 
n = 4 to 1 transition lines were clearly detected, which is a 
firm evidence for charge-exchange emission. We consider 
that the emission is due to the charge-exchange interac- 
tion of solar-wind heavy ions with neutrals in the mag- 
netosheath at 2-8 Earth radii (i?©). In this paper, we 
will concentrate on the X-ray spectra and the emission 
processes, and their implications for cosmic X-ray obser- 
vations. The geophysical implications of the results will 
be discussed in a separate paper. 

Errors quoted in the text and tables are at 90% confi- 
dence single parameter errors and at la confidence level 
in the figures, unless otherwise stated. 

2. Analysis and results 

The NEP region was observed with Suzaku twice during 
the Science Working Group (SWG) observation time. The 
XIS was set to normal clocking mode and the data format 
was either 3 x 3 or 5 x 5. A log of the observations is shown 



in table 1. In this paper, we concentrate on the spectral 
change during a "fiare" detected in the first observation. 
For that purpose, we try to model the stable components, 
and then evaluate the spectral change. We use the data 
from the backside-illuminated CCD (XISl), because of 
its much superior performance below 1 keV (Koyama et 
al. 2006). 

2.1. Data reduction 

We used version 0.6 processed Suzaku data^. In addi- 
tion to the standard data selection criteria: elevation from 
sunlit earth rim > 20°, elevation from dark earth rim > 5°, 
we applied cutoff rigidity (COR) > 8 to clean the XIS data. 
The XIS pulse height data for each X-ray event (3 x 3 or 
5x5 format) were converted to PI (Pulse Invariant) chan- 
nels using the 'xispi' software version 2005-12-26 and CTI 
parameters from 2006-01-25. 

We first created a time history of the X-ray counting 
rate by binning the event data into 256 s time bins. In 
figure 1, we show the 0.3-2 keV counting rate of XISl 
where the non- X-ray (particle- induced) background rate is 
not subtracted. The counting rate shows a clear enhance- 
ment in the first ^ 4 x 10** s. The particle background 
rate of the XIS is known to vary because the cosmic ray 
flux changes as a function of the spacecraft position over 
the Earth. The background rate is well reproduced as a 
consistent function of the local cutoff rigidity. From the 
XIS data during intervals when the telescope is pointed at 
the dark side of the Earth, we found that the 0.3-2 keV 
XISl non- X-ray counting rate varies only from 0.03 cts 
to 0.07 ctss^^ when the cutoff rigidity varies from 15 to 
6 GV. Thus the observed enhancement cannot be particle 
background variation. 

The enhancement lasted for 10 hours. Within the 
enhancement there are shorter time variations. For exam- 
ple, there are sharp spikes just before and after the highest 
peak at ^ 2.2 x lO'' s. A time scale of these variations is 
as short as ^ 10 minutes. We defined the "fiare" interval 
to be to 4 X 10^ s and the "stable" interval as 4 x 10"^ 
to 16.5 X 10"* s as shown in figure 1, then created X-ray 
images for "fiare" and "stable" periods separately. The 
images in figure 2 show that the enhancement is not due 
to changes in any discrete sources. 

We subtracted the particle background from both "sta- 

^ Version processing is an internal processing applied to the 
Suzaku data obtained during the SWG phase, for the purpose of 
establishing the detector calibration as quickly as possible. Some 
processes that are not critical for most of the initial calibration 
and scientific studies, e.g., aspect correction, fine tuning of the 
event time tagging of the XIS data, are skipped in version 
processing, and hence, quality of the products is limited in these 
directions, compared with the official data supplied to the guest 
observers. As of 2006 July, version 0.7 is the latest, where the 
absolute energy scale accuracy of ±0.2 eV at the iron Ko energy 
and ±5 eV below 1 keV is achieved for the XIS data (Koyama 
et al. 2006). In this paper, we used version 0.6 data where the 
energy scale of the XIS data are less accurate (~ 10 eV below 
1 keV) than that of version 0.7, because the empirical model of 
the contamination distribution was obtained based on version 
0.6 data. Instead, we adjusted the scale and the offset of the 
response matrix, as shown in section 2.1. 
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target coordinates 




{a, 


S) = (272.8000, 66.0000) 




observation ID 




100018010 


500026010 




observation period 


2005 


Sep. 2 14:30-Sep 


. 4 15:00 2006 Feb. 10 5:50-Feb 


12. 2:00 


net exposure time 




109.8 ks 


88.6 ks 






(3 X 


3: 95.7 ks, 5x5: 


14.1 ks) (3 X 3: 71.5 ks, 5x5: 


17.1 ks) 



flare 



stable 



r mm 





5x10' 10° 1.5x10= 

sec from 14:54:20, Sep 2, 2005 (UT) 



Fig. 1. XISl counting rate in the 0.3-2 keV energy band (top 
panel), solar-wind proton flux (seeond), flux (third), and 
O"*"^ flux (bottom) as a funetion of time. The Suzaku XIS 
eounting rate is shown in 256 s bins. Particle background 
counts are not subtracted. The solar-wind proton flux was 
calculated using level 2 ACE SWEPAM data. Each bin of 
the ACE data was flrst shifted in time to correct for the 
travel time of the solar wind from ACE to the Earth (typ- 
ically ~ 2700 s), then rebinned into 256 s bins. The ion fluxes 
{C+^, 0+'') were calculated from level 2 ACE SWICS data. 
Only good data with quality flag were used. See also text 
in section 3. 
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Fig. 2. XISl images in the 0.3-2 keV band for "flare" and 
"stable" periods. 
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ble" and "flare" spectra using the dark earth data with 
the same distribution of cutoff rigidities. The count- 
ing rates of the background subtracted spectra above 
10 keV, where the data are dominated by non X-ray back- 
ground events, are (2.9±0.5) x 10~^ ctss~^ for "Hare" and 
(1.2 ±0.3) X 10~2 ctss"^ for "stable", where the error is 
la statistical error. These correspond to 6% and 2% of 
the counting rates of the dark earth data, respectively. 
Note that, in the 2-5 keV band, change of the spectral 
model due to this background difference is much smaller 
than the statistical error, and hence, they do not affect 
the results in the later sections. 

Since we suppose the diffuse X-ray emission to be ap- 
proximately uniform over the field of view of the XIS, we 
created an X-ray telescope (XRT) response function, or 
more specifically, an ancillary response file (ARF) used 
in the XSPEC spectral fit software (Arnaud 1996), for a 
flat field, using the XIS ARF builder software xissimarf- 
gen (Ishisaki et al. 2006). It is known that contamination 
has been accumulating on the optical blocking filters of the 
XIS sensors since the detector doors were opened following 
the launch, and that it accumulates much more quickly at 
the center of the field of view than at the outside (Koyama 
et al. 2006). We used the contaminant thickness and 
radial distribution functions version 2006-05-28 when we 
built an ARF with xissimarfgen. At the early time of this 
NEP observation, the contaminant column densities were 
only 4.1 X 10^^ carbon atoms cm~^ and 0.7 x 10^^ oxy- 
gen atoms cm~^ at the center for XISl. This reduces the 
efficiency by about 12% at the energy of Cvi (367 eV), 
6% at O VII (561 eV), 4% at O Vlll (653 eV), and < 2% at 
NelX (905 eV) and higher energies. Systematic errors in 
the contaminant thickness are estimated to be about 10%. 
The transmission uncertainty due to this systematic error 
is only 1% for Cvi, and less for lines at higher energies; 
hence it is negligible compared with other errors. For the 
XIS response function, wc used the ae_xil_20060213c.rmf 
file supplied by the XIS team, with energy scale correc- 
tions of slope 0.9948 and offset —0.0035 keV as determined 
through the iterative analysis described in the next sec- 
tion. 

2.2. Spectral fit of "stable" spectrum 

We then performed spectral fits to the "stable" spec- 
trum. We first restricted the fitting energy range to 
2-5 keV. In this range the emission is dominated by 
the Cosmic X-ray Background (CXB), which is largely 
emission from unresolved AGNs and can be represented 
by a power law function absorbed by neutral material 
along the line of sight through our Galaxy. Wc thus 
fitted the spectrum with a power-law function with ab- 
sorption by a neutral medium with solar abundances 
(Anders, Grevesse 1989) and fixed the absorbing col- 
umn density at the total Galactic value in this direction 
A^H = 4.4 x 10^0 cm-2 (Dickey, Lockman 1990)^. The fit 
results are summarized in the second column of table 2. 

^ We used 'nH' tool available at http:/ /heasarc.gsfc. nasa.gov/cgi- 
bin/Tools / w3nh/w3nh.pl. 
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Fig. 3. Spectral fit to the "stable" spectrum in the 0.2-2 keV 
range. Observed spectrum is plotted in the upper panel with 
crosses where the vertical bars correspond to la statistical 
errors. The thick step function is the best-fit model function 
convolved with the X-ray mirror and the detector response 
functions. The dotted lines show contributions of different 
spectral components. The lower panel shows the residuals of 
the fits. 

The photon index is consistent with the nominal CXB 
value (1.40 ±0.05; MarshaU et al. 1980). The normaliza- 
tion of the power-law function is also consistent with pre- 
vious observations; 9-11 photons cm~^ s~^ sr~^ keV^^ at 
1 keV (e.g., Gendreau et al. 1995; Revnivtsev et al. 2005). 

We then fitted the spectrum over the 0.2-2 keV range. 
As shown in figure 3, it is clear that there is additional 
emission below 1 keV and the excess contains emission 
hues such as OVII, OVIII and NelX Ka. Thus, fixing 
the spectral parameters of the CXB component to the 
best fit values of the 2-5 keV fit, we added a thin ther- 
mal emission component using the MEKAL model (Mewe 
et al. 1985; Mewc et al. 1986; Kaastra 1992; Liedahl et 
al. 1995). We fixed the abundance of O at the solar value, 
and set abundances of other elements (C, N, Ne, Fe) free. 
The version 0.6 XIS data products are known to contain 
systematic energy calibration errors of ~ 10 eV amplitude 
below ~ 2 keV. We therefore adjusted the energy scale by 
varying the gain slope and offset as free parameters of the 
fit. We determined the gain and the offset in this fitting, 
and applied the same gain and offset throughout the pa- 
per, including the fitting of the CXB component described 
in the previous paragraph. The results are shown in the 
third column of table 2. The model adequately represent 
the observed spectrum'^. Therefore, we adopt the model 
shown in table 2 as a representative model for the "stable" 
spectrum, in order to evaluate the spectral change during 
the "flare". 

The spectrum could be fit by other thermal models. 

^ There are small positive residuals, especially at around 300 eV 
and 450 eV. We can model these features with additional delta 
functions of variable energy and amplitude. Even if we add 
them, however, normalizations of the delta functions employed 
for the "flare" spectrum (section 2.3) are affected by 10% at 
most. 
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Table 2. Results of spectral fits to the "stable" spectrum 
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Parameter 


2-5 keV 


0.2-2 keV 


Power-law component 


TVh [10^^ cm-^] 


0.044 (fixed) 


0.044 (fixed) 


photon index F 


1.33 ±0.23 


1.33 (fixed) 


normalization* 


10.4+3-1 


10.4 (fixed) 



Thin-thermal component (VMEKAL) 



kT [keV] 

C abundance (solar) 

N abundance (solar) 

O abundance (solar) 

Ne abundance (solar) 

Fe abundance (solar) 

normalization^ 
gain slope 0.9948 (fixed) 0.9948 

gain oflFset -0.0035 (fixed) -0.0035 

xVdegrees of freedom 38.39/38 280.15/228 " 

* In units of photons cm^"^ s^^ sr^^ keV^^ at 1 keV. 

t (47r)~^Z?^^(l + 2)~^10~i'' J' ncnudV per steradian, -where Dp^ is the 
angular size distance to the source (cm), and ric, njj are the electron and 
hydrogen densities (cm"-^), respectively. 



-1-0. oua 

0.002 
+0.40 



0.177 

2 14-HO.33 

^■'-^-0.31 

1.0 (fixed) 

'-0.59 

1 42+0-20 

1^-^^-0.22 
iO.O0_Q fig 



When we adopt the present single temperature model with 
varying abundances, the temperature is determined pri- 
marily by the OVII to OVIII Ka emission line intensity 
ratio, and the abundances are determined by the intensi- 
ties of C VI Ka, NVI Ka, FeXVII-L and NelX Ka. If we 
employ a multi-temperature thermal model, it may not be 
necessary to vary abundances because line intensity ratios 
can be adjusted by choosing temperatures. The most im- 
portant result here is that the excess emission above the 
CXB below 1 keV can be represented by the emission lines 
of a thermal model. We tried to include additional con- 
tinuum emission represented by a power-law function or 
thermal bremsstrahlung model. However, there was no 
improvement in . 

2.3. Spectral fit of ''flare" spectrum 

We first compared the "flare" spectrum and the best 
fit model for the "stable" spectrum in the 2-5 keV range. 
We found that there is an excess hard emission (x^/d.o.f 
= 46.21/16) and added an additional power law compo- 
nent, as shown in figure 4. Then the "flare" 0.2-2 keV 
spectrum was compared with a model consisting of the 
best fit "stable" model and the additional power law com- 
ponent as shown in figure 5. The residuals at the bottom 
of the figure show line-like structures, so we have added 
nine emission lines as indicated by the arrows. All lines 
were modeled by delta functions of variable energy and 
amplitude. We show the results in the third column of 
table 3 and figure 6. The "fiare" spectrum is well repre- 
sented by the model"*. Therefore, the enhancement of the 
X-ray intensity during the "flare" can be explained by an 

* The least significant line is that at 882 eV. By adding this line, 
X^/d.o.f -was improved from 188.61/116 to 161.85/114. This line 
is significant at a greater than 99.98% confidence level based on 
the F-test. 




Additional power-law component 
%, required for the "flare" spectrum 

^ 2 3 4 5 

Energy (keV) 



Fig. 4. Comparison of the "flare" spectrum with the CXB 
model fit to the "stable" spectrum at 2—5 keV. An addi- 
tional power-law component was required in this range for 
the "flare" spectrum. Parameters are given in the second col- 
umn of table 3. 

increase in these emission lines and the hard power-law 
emission. 

In the fourth column of table 3, we show probable iden- 
tifications of lines. The lowest energy line below the car- 
bon edge at 269 ± 4 eV is considered to be a sum of multi- 
ple L emission lines. The line at 455^^3 eV is most likely 
the n = 4 to 1 transition (Ly7) of C VI at 459 eV, since 
the line energy is consistent within the statistical error and 
there are no other likely emission lines at this energy. If we 
introduced two lines at 436 eV (C VI Ly/3; n = 3 to 1) and 
459 eV (C VI Ly7) instead of one line of variable energy, 
the normalizations of these lines became 1.93+n nQ and 
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Table 3. Results of spectral fits to the "flare" spectrum. Parameter values of components added to the "stable" spectral model. 



Parameter 


2-5 keV 


0.2-2 keV 


Line identification* 


Additional Power-law co 


mponent without absorption 




photon index P 


0.0±1.7 


0.0 (fixed) 




normalization^^ 




0.74 (fixed) 




Additional narrow Gaussian lines 


1 center energy [eV] 


— 


269 ±4 


C band lines 


normalization^ 




7.7+1-? 




2 center energy [eV] 
normalization^ 


— 



3571^ 
7.3t?-^ 


CVI 2p to Is (367 eV) 


3 center energy [eV] 
normalization-!- 


— 


455+^3 


CVI 4p to Is (459 eV) 


4 center energy [eV] 
normalization-!- 


— 


558+^ 
5 1+1-1 

"3-1-1.0 


OVII (561 eV) 


5 center energy [eV] 




649+^ 


OVIII 2p to Is (653 eV) 


normalization-!- 




r 02+0.58 




6 center energy [eV] 
normalization^ 




796+1" 

l-D'-0.34 


Fexvil,xvill-L + Ovill 3p to Is (774 eV)? 


7 center energy [eV] 
normalization-!- 




882+}^ 


Fexvil,xvill-L + Neix (905 eV) 
+ OVIII 6p to Is (847 eV)? 


8 center energy [eV] 




1022+^1 


Nex (1022 eV) 


normalization-!- 




1 n4+"-2o 

1-^^-0.29 




9 center energy [eV] 
normalization-!- 




1356+20 
73+0-19 

^- ' "J-0.20 


Mgxi (1329 eV) 


/d.o.f 


17.65/14 


161.85/114 





* Line energies at the rest frame are taken from Kharchenko et al. (2003), Krasnopolsky (2004). Energies of the forbidden 

line are shown for Ovil, Neix, and Mgxi Ka, because the forbidden line becomes much stronger at the charge-exchange 

emission (e.g., Kharchenko et al. 2003). 

t In units of photons cm~-^ s~l sr~l keV~^ at 1 keV. 

■f In units of photons cm^^ s^i sr~i. 
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Fig. 5. Comparison of the "flare" spectrum with the best flt 
model for the "stable" spectrum plus additional power-law 
component during the "flare" in the 0.2-2 keV energy range. 
Vertical arrows indicate line-like structures in the residuals. 



Fig. 6. "Flare" spectrum and best-fit model spectrum. The 
model spectrum is a sum of the best fit model for the "stable" 
spectrum plus additional power-law component during the 
"flare" shown with a thin solid curve (same as that shown in 
figure 5), and nine emission lines shown with dotted curves. 
Best-fit parameters of the nine emission lines are summarized 
in table 3. 
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1.68^Q Y^ photons cm^^ s^^ sr^^, respectively. Therefore, 
C VI Ly/? could have a comparable contribution. In ci- 
ther case, wc definitely need the C VI Ly7 line. Dcnnerl et 
al. (2003) attributed a weak peak structure found in the 
XMM-Newton spectrum of comet C/2000 to a sum of C VI 
Ly/3 and Ly7, and C VI Ly7 due to charge-exchange emis- 
sion between the highly ionized solar wind and exospheric 
or interplanetary neutrals during an XMM-Newton obser- 
vation of the Hubble Deep Field-North was reported by 
Snowden, Collier, Kuntz (2004). This NEP observation, 
however, seems to be the clearest detection so far. We 
also have detected OVII to MgXI lines. The fines at 796 
and 882 eV are likely to represent complex structures due 
to Fe-L and other lines. 

3. Discussion 

The short (~ 10 minutes) time scale variations observed 
during the enhancement of X-ray intensity imply that the 
size of the emission region is no larger than 10 light min- 
utes. On the other hand, the apparent size of the emis- 
sion region must be equal to or larger than the XIS field 
of view (18'). These require the emitter of the X-ray en- 
hancement to be within a distance of 10 light minutes/18', 
or ~ 10~^ pc. Because the enhanced X-ray emission con- 
sists of emission lines from C VI to MgXI, this requires an 
ion source within 10~^ pc. There is only one ion source 
in this distance range. That is the Sun. 

The Sun may produce X-ray emission lines in our ob- 
servations in two possible ways: scattering of solar X- 
rays by the Earth's atmosphere, and solar-wind charge ex- 
change. In the former case, the X-ray intensity is propor- 
tional to the solar X-ray intensity multiplied by the sunlit 
atmospheric column density. The solar X-ray intensity 
is continuously monitored by the GOES (Geostationary 
Operational Environmental Satellites)^, but data show no 
correlation with the enhancement. Moreover, using the 
MSIS atmosphere model (Hedin 1991), we found that the 
column density of sunlit atmosphere varied by many or- 
ders of magnitude (10^) during the observations, but no 
correlation was found with the observed X-ray intensity. 
Thus scattering of solar X-rays can be excluded. 

In figure 1, we show the proton fiux observed by the 
ACE (Advanced Composition Explorer)^ together with 
the Suzaku X-ray counting rate. The ACE data were 
shifted in time to account for propagation time from ACE 
to the Earth. Clearly the proton fiux was enhanced dur- 
ing the X-ray "flare" . This is consistent with solar-wind 
charge-exchange model. 

Charge-exchange X-ray emission also is strongly sup- 
ported by the detection of the Cvi ti = 4 to 1 transi- 
tion line (Ly7). In CIE (coUisional ionization equilibrium) 
thermal emission, the Ly/3 and Ly7 lines of C VI are sup- 
pressed relative to Lya by the Boltzmann factor in the 
distribution of exciting electrons. In charge exchange be- 
tween C VII and Hi, the electron is deposited primarily 

^ The data available at http:/ /www. ngdc.noaa.gov/stp/GOES/goc 
® The data available at http://www.srl.caltech.edu/ACE/ASC/. 
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Fig. 7. 0.3-2 keV X-ray intensity (upper panel) and geocen- 
tric distance of the point whose geomagnetic field becomes 
open to the space for the first time along the line of sight 
from the spacecraft position in units of the earth radius, as 
a function of time during the observation (lower panel). See 
also the schematic view shown in figure 8. 
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Fig. 8. Schematic view of the magnetosphere and the defini- 
tion of rmp used in figure 7. 

in the n = 4 level (Krasnopolsky et al. (2004) and ref- 
erences therein), and the X-ray lines arc produced when 
it cascades to the n = 1 level, guided only by branch- 
ing ratios. In high energy collisions, angular momentum 
states are tend to be populated statistically by weight of 
the state's degeneracy, and electrons are primarily cap- 
tured into maximal /, where n can change only by 1 unit 
at a time during the cascade (Beiersdorfer et al. 2001), 
again resulting in relatively weak Ly/3 and Ly7. Behind 
the Earth's bow shock, however, the solar wind velocity is 
reduced, which should result in recombination into low l- 
orbitals and strong Ly/3 and Ly7. In the "fiare" spectrum, 
both Cvi Lya and Ly7 lines were detected and this is a 
strong evidence for charge exchange. In addition, O VIII 
n = 3 to 1 and rt = 6 to 1 transition lines are also enhanced 
in^jthc comet emission model (Kharchenko et al. 2003). 
However, we cannot separate those lines from FeXVII-L 
emission lines with the present energy resolution. 
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Although X-ray intensity and the proton flux are cor- 
related on a time scale of ^ 10 hours, they do not show 
much correlation on short time scales. We consider that 
the short term X-ray intensity variation is at least partly 
arising from orbital motion of the spacecraft. In figure 7, 
we show the geocentric distance of the point whose geo- 
magnetic field becomes open to the space for the first time 
along the line of sight from the spacecraft position, i.e., 
the point where the line of sight encounters the magne- 
tosheath (see also the schematic view shown in figure 8 for 
the definition). We evaluated the end point of the mag- 
netic field using the software GEOPACK-2005 and T96 
magnetic field model (Tsyganenko, Sitnov (2005) and ref- 
erences therein) We obtained the solar- wind parameters 
required to perform the calculation from the CDAWeb 
(Coordinated Data Analysis Web)^. We find that the line 
of sight during the present observation was rather spe- 
cial in the sense that it goes through the north magnetic 
pole region where charged particles of the magnetosheath 
can penetrate down to 2-8i?0 moving along open field 
lines. The short-term X-ray intensity variations during 
the time intervals shown by boxes in figure 7 indicate 
anti-correlation with the distance to the magnetosheath. 
This indicates that the charge exchange of solar-wind 
heavy ions is taking place at 2-8i?0 where the neutral 
matter density is high. Robertson et al. (2006) recently 
studied theoretically solar- wind charge- exchange emission 
from the magnetosheath. Implications of the present re- 
sults on the solar-wind ion composition and the Earth's 
magnetosheath including comparisons with the theoreti- 
cal model will be reported in a separate paper. 

Finally, since solar-wind charge- exchange emission can 
become a difficult foreground in the study of soft diffuse 
sources, we summarize a procedure to examine possible 
contamination in the Suzaku spectra below 1 keV by 
solar activity. 

1. check the light curve, if time variation is not ex- 
pected for the object, 

2. check the solar X-ray intensity and the column den- 
sity of sunlit atmosphere along the line of sight, 

3. check the solar wind proton flux, 

4. check the radius of the magnetosheath on the line of 
sight. 

We would like to thank Prof. K. Maezawa, and 
Prof. I. Shinohara for their help in calculations of the 
Earth's magnetosheath and also for valuable discussions. 
Thanks are also due to Prof. T. Mukai, Prof. A. Matsuoka, 
and Prof. H. Hayakawa for discussions on the relation be- 
tween Suzaku data and solar wind. We are grateful to the 
referee. Dr. A. Dalgarno, for useful comments to improve 
this paper. 

^ The software package available at 

http: / /modolwcb.gsfc.nasa.gov/magnctos/data- 
bascd/modcling.html. 

® The data available at http:/ /cdaweb. gsfc.nasa.gov/cdawcb/sp_pl] 
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